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Abstract 

\ Radiative transfer in a relativistic plane-parallel flow, e.g., an accretion disk wind, is examined in the 

fully special relativistic treatment. Under the assumption of a constant flow speed, for the relativistically 
moving atmosphere we analytically obtain generalized Milne-Eddington solutions of radiative moment 
equations; the radiation energy density, the radiative flux, and the radiation pressure. In the static limit 
these solutions reduce to the traditional Milne-Eddington ones for the plane-parallel static atmosphere, 
whereas the source function nearly becomes constant as the flow speed increases. Using the analytical 
00 ' solutions, we analytically integrate the relativistic transfer equation to obtain the specific intensity. This 

specific intensity also reduces to the Milne-Eddinton case in the static limit, while the emergent intensity 
is strongly enhanced toward the fiow direction due to the Doppler and aberration effects as the flow speed 
pq ' \ increases (relativistic peaking) . 
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1. Introduction 

^ ' Radiative transfer problems in astrophysics have been very widely and extensively studied since 1920's (Milne 1921; 
Eddington 1926; Kosirev 1934; Chandrasekhar 1934). There are many excellent monographs on radiative transfer and 
radiation hydrodynamics (e.g., Chandrasekhar 1960; Mihalas 1970; Rybicki, Lightman 1979; Mihalas, Mihalas 1984; 
\ Shu 1991; Peraiah 2002; Castor 2004). In the early days the researchers were focused their attention on radiative 
. transfer in a static atmosphere of plane-parallel and spherical cases, relating to stellar atmospheres and nebulae. Later 
they examined radiative transfer in extended and expanding atmospheres, relating to stellar winds, nova explosions, 
and so on. 

Todays we know there are many energetic outflows in the universe, which are blown off from the central luminous 
sources with mildly to highly relativistic speeds. For example, accretion disk winds can emanate from cataclysmic 
variables (CVs), supersoft X-ray sources (SSXSs), X-ray binaries (XBs), and active galactic nuclei (AGNs); relativistic 
\ jets are ejected from microquasars (^QSOs) and quasars (QSOs); relativistic outflows occur in supernovae (SNe), 
hypernovae (HNe), and gamma-ray bursts (GRBs). 

In such circumstances, radiative transfer problems in relativistically moving media become more and more important. 
However, there are little theoretical and/or analytical studies. In particular, radiative transfer in the accretion disk 
wind has not been well considered. For example, transformation properties of disk radiation fields in the proper 
frame of a relativistic jet were examined by, e.g., Dermer and Schlickeiser (2002). In these earlier works, however, 
the radiation fields are set to be external sources, and the radiation transfer was not considered. Recently, radiative 
transfer in a moving disk atmosphere was firstly investigated in the subrelativistic regime (Fukue 2005a, 2006a), and 
in the relativistic regime (Fukue 2005b, 2006b; Fukue, Akizuki 2006b). In these studies, however, only the radiative 
moments were obtained under the moment formalism, and the specific intensity was not solved. Hence, the specific 
intensity in a plane-parallel moving atmosphere was obtained in the subrelativistic regime (Fukue 2007), and in the 
relativistic regime (Fukue 2008a). 

In the previous study (Fukue 2008a), under the assumptions that both the fiow speed and the source function are 
constant, we were able to obtain the analytical solutions. Since we assumed that the source function is constant, 
however, the solutions obtained cannot be compared with the traditional solutions such as Milne-Eddington ones for 
a static atmosphere. 

Thus, in this paper we loosen the assumptions and extend the previous work. Namely, we examine radiative transfer 
in the relativistic plane-parallel fiow, such as an accretion disk wind, which is assumed to blow off from the luminous 
disk in the vertical direction (plane-parallel approximation) , at a constant speed. 

In the next section we describe the basic equations. In sections 3 and 4, we show analytical solutions of moment 
equations, and of transfer equation, respectively. We then briefly discuss the consistency of the present analytical 
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solutions in section 5. The final section is devoted to concluding remarks. 

2. Relativistic Radiative Transfer Equation and Moment Equations 

Let us suppose a plane-parallel moving atmosphere, e.g., a luminous flat disk wind. The radiation energy is 
transported in the vertical direction, and the gas, itself, also moves in the vertical direction as a plane-parallel outflow 
due to the action of radiation pressure. For simplicity, in the present paper, the radiation field is considered to be 
sufficiently intense that both the gravitational field of, e.g., the central object and the gas pressure can be ignored. 
We also assume the gray approximation, where the opacities do not depend on the frequency. As for the order of the 
flow velocity v, we consider the fully special relativistic regime. 

The radiative transfer equations are given in several literatures (Chandrasckhar 1960; Mihalas 1970; Rybicki, 
Lightman 1979; Mihalas, Mihalas 1984; Shu 1991; Kato et al. 1998, 2008; Mihalas, Aucr 2001; Peraiah 2002; Castor 
2004). The basic equations for relativistic radiation hydrodynamics are given in, e.g., the appendix E of Kato ct al. 
(2008) in general and vertical forms, where the quantities are expressed in the inertial frame. The relativistic radiation 
hydrodynamic equations in the laboratory frame are also derived and discussed in the seminal paper by Mihalas and 
Auer (2001) in some details. 

In a general form the radiative transfer equation in the inertial (fixed) frame is expressed as 
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Here, v is the flow velocity, c is the speed of light, and 7 (= 1/^1 — jc?') is the Lorcntz factor. In the left-hand side 
the frequency-integrated specific intensity J and the direction cosine I are quantities measured in the inertial (fixed) 
frame. In the right-hand side, the mass density p, the frequency- integrated mass emissivity jo, the frequency- integrated 
mass absorption coefficient , and the frequency-integrated mass scattering coefficient k^q'^ are quantities measured 
in the comoving (fluid) frame, whereas the frequency-integrated radiation energy density i?, the frequency-integrated 
radiative flux F , and the frequency-integrated radiation stress tensor P'-' are quantities measured in the inertial (fixed) 
frame. 

In the plane-parallel geometry with the vertical axis z and the direction cosine /i (= cos0), the transfer equation is 
expressed as 
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where (3 (= v/c) is the normalized vertical speed, and F and P are the vertical component of the radiative flux and 
the radiation stress tensor measured in the inertial frame, respectively. 

For matter, the continuity equation, the equation of motion, and the energy equation become, respectively. 
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where u (= 7/3) is the vertical four velocity, J the mass-loss rate per unit area, ip the gravitational potential, e the 
internal energy per unit proper volume, p the gas pressure, and the internal heating. In the energy equation (5) 
the advection terms in the left-hand side are dropped under the present cold approximation (radiative equilibrium). 
We do not use the equation of motion since we will assume that the flow speed is constant. 

For radiation, the zeroth and first moment equations, and the closure relation become, respectively, 
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where /(t, /?) is the variable Eddington factor, which is defined by / = Pco/^E'coj -Sco and Pco being the comoving 
quantities, and generally depends on the velocity and its gradient as well as the optical depth (Fukue 2006b, 2008b). 
In the present approximation, however, / is automatically determined to be 1/3 as will be shown later. 
Eliminating jo using the energy equation (5), the transfer equation (2) becomes 
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Introducing the optical depth defined by 
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the transfer equation (9) finally becomes 
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is the scattering albedo, and the internal heating is dropped. 

Similarly, moment equations (6) and (7), with the help of continuity equation (3) and the closure relation (8), 
become (cf. Fukue 2005b, 2006b) 
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Here, we dropped the gravitational and pressure forces as well as the internal heating. 
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3. Solutions of Moment Equations 



In the present paper we assume that the flow speed is constant, but the source function is not, and seek the analytical 
solutions of moment equations (13) and (14). We also assume that the Eddington factor docs not depend on the optical 
depth. 
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3.1. General Solutions 
If equation (13) is divided by equation (14), we have 
dF 



(15) 



which is integrated to give 

F = /3cP + Fi(/3), (16) 

where Fi (P) is an arbitrary function of (3. 

Inserting this relation (16) into equations (13) and (14), we can easily integrate equations (13) and (14) to give the 
following forms: 
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where Cf(/3) and Cp(/3) are arbitrary functions of (3. 

In order to both solutions (17) and (18) are consistent with each other, Cp = Cp/3 and we have 
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and furthermore, from the closure relation (8), we also have 
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With the help of arbitrariness in Cp, we can express general solutions (17), (19), and (20) as 
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where Ci(/3) = 7^(/ — /3^)FiCp is an arbitrary function of /?. 

These equations (21)-(23) are general solutions of rclativistic moment equations for plane-parallel vertical flows at 
a constant speed. In the limit of /3 ^ 0, these solutions approach the Milne-Eddington ones, as long as 1 — Cp ~ /3 for 
small p. 

3.2. Special Solutions 

Next, we derive more special solutions under the appropriate boundary conditions. For the present plane-parallel 
flow at a constant flow speed, the boundary conditions should be imposed at the flow top, where the optical depth is 
zero. 

At the flow top of a moving photosphere at a rclativistic speed, the usual boundary conditions for a static atmosphere 
is inadequate, as already pointed out in Fukue (2005b). Namely, the radiation field just above the flow top changes 
when the gas itself does move upward, since the direction and intensity of radiation change due to the relativistic 
aberration and Doppler effect (cf. Kato et al. 1998, 2008; Fukue 2000). If a flat infinite plane with surface intensity Is 
in the comoving frame is not static, but moving upward at a speed Vs (= C/Sg, and the corresponding Lorentz factor is 
7s), where the subscript s denotes the values at the surface, then, just above the surface, the radiation energy density 
Eg, the radiative flux _Fs, and the radiation pressure Ps measured in the inertial frame become, respectively. 
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It should be noted that these boundary conditions imply that the Eddington factor / defined in the comoving frame 
is 1/3. 

Using these boundary conditions (24)-(26) to general solusions (21)-(23), we can determine the arbitrary functions 
Fi and Ci as 



Fi=^/,(l + 2/3), 



(27) 
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and the Eddington factor / is automatically fixed as 1/3. 

Finally, the solutions of moment equations for relativistic plane-parallel flows at a constant speed are explicitly 
written as 
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where / = 1/3. 

In the limit of small /3, these solutions approach 
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which are just the Milne-Eddington solutions for a static plane-parallel atmosphere. Hence, the present solutions 
(29)-(31) are the Milne-Eddington ones extended to the case for the relativistic plane-parallel flow. 

By the Lorentz transformation, it is easy to obtain the solutions in the comoving frame. Namely, the solutions of 
moment equations in the comoving frame are explicitly written as 
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where / = 1/3. Hence, as already stated, in the comoving frame Pco/Eco — f — 1/3, exactly. 

In figure 1, Milne-Eddington solutions of moment equations for relativistic plane-parallel flows are shown as a 
function of the optical depth r: (a) Radiation energy density E, (b) radiative flux F, and (c) radiation pressure P, 
where the radiation quantities are normalized by tt/s. The solid curves represent the quantities in the inertia! frame, 
whereas the dashed ones mean those in the comoving frame. The flow speed (3 are set to be 0, 0.1, 0.2, 0.4. and 0.577. 
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Fig. 1. Milnc-Eddington solutions for relativistic plane-parallel flows: (a) Normalized radiation energy density, (b) normalized 
radiative flux, and (c) normalized radiation pressure. The solid curves represent the quantities in the inertial frame, whereas the 
dashed ones mean those in the comoving frame. The values of f3 are 0, 0.1, 0.2, 0.4, and 0.577. 

In figure la, the radiation energy density E normafized by ttIs/c is plotted for various f3. As was already stated, 
and as is easily seen in figure la, when the flow speed is small, the solution reduces to the usual Milne-Eddington one 
for a static plane-parallel atmosphere; cE/{ttIs) ^ 2 + 3t. As the flow speed increases, the radiation energy density 
becomes constant. This means that the radiation energy itself is advected along with the flow, and it is constant 
without the internal heating. Furthermore, comparing the radiation energy density in the inertial frame (solid curves) 
and those in the comoving frame (dashed ones), the radiation energy density in the inertial frame is enhanced due to 
the relativistic motion. 

In figure lb, the radiative flux F normalized by Tr/g is plotted for various /3. When the flow speed is small, the 
solution reduces to the usual Milne-Eddington one for a static plane-parallel atmosphere; F/{ttIs) ^ 1. As the flow 
speed increases, the radiative flux in the inertial frame (solid curves) is enhanced more and more due to the relativistic 
effect. The radiative flux in the comoving frame (dashed ones), on the other hand, reduces exponentially. This means 
that in the rclativistically moving plane-parallel flows without the internal heating the radiative flux in the comoving 
frame vanishes since the radiation energy (and pressure) is advected with the flow and there is no gradient in the 
radiation pressure. 

In figure Ic, the radiation pressure P normalized by Tr/g/c is plotted for various (3. When the flow speed is small, 
the solution reduces to the usual Milne-Eddington one for a static plane-parallel atmosphere; cP/{'kIs) ^ 2/3 + t. 
As the flow speed increases, the radiation pressure becomes constant. Furthermore, comparing the radiation energy 
pressure in the inertial frame (solid curves) and those in the comoving frame (dashed ones), the radiation pressure in 
the inertial frame is enhanced due to the relativistic motion. 

It should be noted that in the present solution the Eddington factor / is set to be 1/3, and therefore, the flow speed 
(3 is restricted within 13 < 1/^3-0.5774 (cf. Fukue 2005a). 
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4. Solutions of Transfer Equation 

Since we obtain the analytical forms of radiation moments, and we assume that the flow speed is constant, we can 
analytically integrate the relativistic radiation transfer equation (11). 

By inserting the solutions (29)-(31) into the transfer equation (11), and setting / = 1/3, the transfer equation (11) 
is rewritten as 



(38) 



where / = 1/3. 

It is straightforward to integrate this equation (38) (cf. Fukue, Akizuki 2006b; Fukue 2007). After several partial 
integrations, we obtain both an outward intensity /(t, /i,/3) (/i > 0), 
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where I{tq,^) is the boundary value at the flow base of the optical depth tq. 

In the limit of small (3, this solution reduces to the usual Milne-Eddington one; e.g., for /i > 
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Finally, the emergent intensity /(0,/i,/3) (// > 0) at the flow surface of r = becomes 
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Fig. 2. Normalized specific intensity of Milne-Eddington solutions for relativistic plane-parallel fiows as a function of the optical 
depth T for several values of fi and /3. The values of /3 are (a) 0, (b) 0.1, and (c) 0.3. The values of fi arc (thick solid curves), 
0.2,0.4,0.6,0.8, and 1 (solid curves), and -0.2,-0.4,-0.6,-0.8, and -1 (dashed ones). 



^Is 1 ^"^7'(1-/3m)(/-/^') 
An 74(1 - /3/,)4 



for large tq. (42) 



In figure 2, the specific intensity normalized by Is of Milne-Eddington solutions of transfer equation for relativistic 
plane-parallel flows are shown as a function of the optical depth r for several values of n and /3. The values of P are (a) 
0, (b) 0.1, and (c) 0.3. In each figure, the thick solid curves are for n = 0, the solid curves are for fi = 0.2,0.4,0.6,0.8, 
and 1 from bottom to top, while the dashed ones are for n = —0.2, —0.4, —0.6, —0.8, and —1 from top to bottom. 

As is seen in figure 2, in the case of small (3 (figure 2a), the usual Milne-Eddington solution for a static plane-parallel 
atmosphere is well reproduced. As the flow speed increases, the specific intensity becomes constant, and is enhanced 
(/i > 0) or diminished (^ < 0) due to the relativistic vertical motion. 

In figure 3, the emergent intensity /(0,/i,/?) normalized by Is of Milne-Eddington solutions of transfer equation for 
relativistic plane-parallel flows is plotted as a function of fi for several values of /3. The values of (3 are 0, 0.1, 0.2, 0.3, 
0.4, and 0.5. The dashed line is for the usual Milne-Eddinton solution for the static plane-parallel case. 

In the limit of /3 = 0, the present emergent intensity again reduces to the usual Milne-Eddington solution for a static 
atmosphere. As the flow speed increases, the relativistic peaking effect becomes prominent; that is to say, the emergent 
intensity is strongly enhanced toward the poleward direction, while it is less enhanced in the edgeward direction. 

5. Discussion 

In the previous sections we have obtained the analytical solutions (29)-(31) of relativistic moment equations, and 
the analytical solutions (39) and (40) of relativistic transfer equation, using the solutions (29)-(31). However, there 
is no assurance that these solutions are consistent each other. Hence, we here briefly discuss the consistency of the 
present solutions, and check the accuracy of the present results. 

Using the radiation fields E, F, and P in the inertial frame, we obtained the radiation intensity / in the inertial frame. 
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Fig. 3. Normalized emergent intensity of Milne-Eddington solutions for relativistic plane-parallel flows as a function of fi for several 
values of /3. The values of (3 are 0, 0.1, 0.2, 0.3, 0.4, and 0.5. The dashed line is for the usual Milne-Eddinton solution for the static 
plane-parallel case. 

Once the radiation intensity /(T,/i,/?) is obtained, we can calculate the radiation energy density Enum (= 27r J Idfi/c), 
the radiative flux Fnum (= Stt J I^dfi), and the radiation pressure Pnum (= 27r / Ifi^dfi/c) by the definition. If the 
quantities obtained from the moment equations and those obtained from the intensity coincide each other, the solutions 
are consistent, and vice versa. 

In figure 4 we show the analytical solutions of moment equations (dashed curves) and the numerical solutions 
obtained by the integration of the analytical solutions of transfer equation (solid curves) for various values of /3: (a) 
Radiation energy density E, (b) radiative flux F, and (c) radiation pressure P, where the radiation quantities are 
normalized by Tr/g. The flow speed /3 are set to be 0, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.577. 

As is seen from figure 4, in the case of small velocity, both solutions are almost coincident. As the flow speed 
increases, the numerical solutions obtained from the intensity become a little bit larger than the analytical solutions. 
However, the discrepancy is not so large, and within ten parcents at most. 

We can also calculate the Eddington factor / by the numerical integration of the analytical solution /(t, /i,/3). In 
figure 5 the values of 3/ are shown for various values of /3 (= 0, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.577). 

In the limit of /3 = 0, i.e., in the case of the static atmosphere, the Eddington factor slightly becomes larger than 
1/3 for the small optical depth. This is just the limb-darkening effect (or the peaking effect) originated from thea 
anisotropic properties of the radiation fields near to the surface. Indeed, it is easy to derive that the Eddington factor 
of the traditional static atmosphere becomes (17/14)(l/3) at r = (see, e.g., Fukue, Akizuki 2006a). 

As the flow speed increases, however, the Eddington factor near to the surface becomes 1/3, and then further 
decreases to be smaller than 1/3. This can be interpreted by the relativistic effect; the relativistic aberration would 
cancell or dominate the limb-darkening effect in the region of small optical depths. Anyway, except for the regions of 
T < 1, the Eddington factor is almost 1/3 even if the flow speed is large. 

6. Concluding Remarks 

In this paper wc have examined the radiative transfer problem in the relativistic plane-parallel flows, e.g., accretion 
disk winds, under the assumption of a constant flow speed in the fully special relativistic regime. We have found the 
analytical solutions (29)-(31) of the relativistic moment equations, and the analytical solutions (39) and (40) of the 
relativistic transfer equation, using solutions (29)-(31). These analytical solutions are the generalized Milne-Eddington 
solutions for the relativistically-moving plane-parallel flows. 

In the subrelativistic regimes many researchers have examined the radiative transfer and radiation hydrodyanmics, 
and the problem of radiation hydrodynamics were often numerically solved (e.g., Ohsuga et al. 2005; Ohsuga 2006; 
see also Mihalas, Mihalas 1984). It seems, however, that the radiative transfer problem in the fully special relativistic 
regime is not yet well understood. Actually, we found new analytical solutions, and there would exist another solutions. 
The analytical approach in the present work as well as the numerical approach would help us to understand the 
relativistic radiative transfer. 

In the present study we assumed that the total optical depth is sufficiently large and the flow speed is constant. If 
the total optical depth is finite, and the boundary condition at the flow base is taken account of, the radiation fields 
and the resultant Eddington factor would be changed (cf. Fukue, Akizuki 2006a, b; Koizumi, Umcmura 2007). Or, if 
the velocity is not constant, but there is a velocity gradient, then the Eddington factor is proved to be less than 1/3 
deep inside the flow (Fukue 2008b). Such cases should be examined in the future. 
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Fig. 4. Radiative quantities obtained from moment equations and those obtained from the radiative intensity: (a) Radiation energy 
density, (b) radiative flux, and (e) radiation pressure. The dashed curves are the analytical solutions of moment equations and the 
solid curves are the numerical solutions obtained by the integration of the analytical solutions of transfer equation. The values of /3 
are 0, 0.1, 0.2, 0.3, 0.4, 0.5, and 0.577. 




Fig. 5. Eddington factor obtained from the radiative intensity for various values of the flow speed. The values of f3 are 0, 0.1, 0.2, 
0.3, 0.4, 0.5, and 0.577. 
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We found that the emergent intensity depends on the flow speed as well as the direction cosine, and exhibits 
a rclativistic peaking effect. As a result, a wind luminosity would be overestimated by a pole-on observer and 
underestimated by an edge-on observer, when we observe an accretion disk wind (cf. Sumitomo et al. 2007). 

We have also assumed the frequency- independent opacities (gray approximation). This is an approximation, and 
can be extended using mean opacities obtained from a frequency-dependent radiative transfer solution (cf. Mihalas, 
Auer 2001). In addition, we have obtained the frequency-integrated intensity /, whereas the frequency-dependent 
intensity I^, should be also examined in the future. 

Furthermore, we have imposed various assumptions, including cold approximation, no heating source, and so on. By 
relaxing these assumptions and integrating the rclativistic transfer equation numerically, we could obtain the emergent 
intensity and spectra more quantitatively. These are also left as future works. 



This work has been supported in part by a Grant-in- Aid for Scientific Research (18540240 J.F.) of the Ministry of 
Education, Culture, Sports, Science and Technology. 
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